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SUMMARY 


Part I contains the results of a preliminary experimental investi- 
gation of a particular design of an underwater ramjet or hydroduct 
powered by compressed air. The hydroduct is a propulsion device in 
which the energy of an expanding gas imparts additional momentum to a 
stream of water through mixing. The hydroduct model had a fineness 
ratio of 5 . 9 , a maximum diameter of 3.2 inches, and a ratio of inlet 
area to frontal area of 0.32. The model was towed at a depth of 1 inch 
at forward speeds between 20 and 60 feet per second for airflow rates 
from 0.1 to 0.3 pound per second. Longitudinal force and pressures at 
the inlet and in the mixing chamber were determined. 

The hydroduct produced a positive thrust -minus -drag force at every 
test speed. The force and pressure coefficients were functions primarily 
of the ratio of weight airflow to free- stream velocity. The maximum 
propulsive efficiency based on the net internal thrust and an isothermal 
expansion of the air was approximately 53 percent at a thrust coefficient 
of 0.10. The performance of the test model may have been influenced by 
choking of the exit flow. 

Part II is a theoretical development of an underwater ramjet using 
air as "fuel." The basic assumption of the theoretical analysis is that 
a mixture of water and air can be treated as a compressible gas. More 
information on the properties of air-water mixtures is required to con- 
firm this assumption or to suggest another approach. A method is sug- 
gested from which a more complete theoretical development, with the 
effects of choking included, may be obtained. An exploratory computa- 
tion, in which this suggested method was used, indicated that the effect 
of choked flow on the thrust coefficient was minor. 


INTRODUCTION 


The basic concept of the hydroduct or underwater ramjet consists 
of the production of thrust by the transfer of the potential energy of 
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a compressed, gas to a flowing liquid through a mixing process. Water 
entering the hydroduct inlet is diffused to a high static pressure and 
mixed with the expanding gas; the mixture is then expanded through the 
duct exit with a resultant increase in total momentum. A theoretical 
analysis of the hydroduct is presented in reierence 1; reference 2 con- 
tains some results of tests of a hydroduct. 

The hydroduct is of potential interest for the propulsion of surface 
craft designed for high speeds because of certain apparent advantages 
over conventional propulsion systems. Because of the limitations of 
planing craft in rough water, any significant increase in the speed of 
practical surface craft must be accomplished with hydrofoils. It has 
been proposed that the hydroduct and the ducts to supply air to it could 
be contained within the hydrofoils and supporting struts. Also, it has 
been conjectured that such a system could have less frontal area, a 
higher speed capability, and lower noise level than supercavitating 
gear-driven propellers; however, an evaluation is not possible at the 
present time because of insufficient hydroduct data. Because no static 
thrust is developed by the hydroduct, the system must include an auxiliary 
propulsion device. 

Because of the apparent advantages of the hydrofoil-hydroduct system 
for propulsion of high-speed surface craft, a preliminary investigation 
of a simple axi symmetric air-powered hydroduct has been made in Langley 
tank no. 1 (this facility has recently been transferred to the David 
Taylor Model Basin) to evaluate the thrust-producing capability of such 
a device. The hydroduct was towed at shallow draft at speeds between 
20 and 60 feet per second with variations in airflow rate from 0.1 to 
0.3 pound per second. Measurements were made of the towing force and 
pressures at the inlet and maximum-area section. The results of these 
tests are described in part I of this report. 

As a result of the preliminary experimental investigation it was 
apparent that a theory was needed to guide future work. Accordingly, 
an analysis was made of the ideal performance of the air-powered hydro- 
duct. This analysis is presented in part II. 


SYMBOLS 


A area, sq ft 

Af frontal area of body, 0.0559 sq ft 


a 


speed of sound, ft/sec 


vo -p" ro h tr 1 
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axial-force coefficient , 


•A 


| PwV» 2 A f 


pressure coefficient, 


P ~ Poo 
■p Pw^°o 




■'T,f 


*A,o 

M 

m 

P 

Pt 


thrust coefficient, ■= 

f p„v» 2 a c 


thrust coefficient based on frontal area, 


axial force in forward direction, lb 


Fiji 


1 IT 2. 

2 Pw^oo A-f 


axial force in forward direction with nose and tail fairings 
on model, lb 

net jet thrust, F ^ - F^ Q , 

Mach number , V/a 

mass flow per unit time, slugs/sec 
static pressure, lb/sq ft 
total pressure, lb/sq ft 


dynamic pressure, , lb/sq ft 


R 

T 

t 

V 


gas constant for air, 53*3 

internal thrust, lb 
water temperature, 505° R 
velocity, ft/sec 


ft -lb 
lb-°F 
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w weight flow, lb/sec 

7 adiabatic exponent (assumed to be 1.4 for air) 

~y average adiabatic exponent 

tj efficiency 

Fp V oo 

experimental efficiency, 

p j 

Rtw a log -a 
■^00 

p density, slugs/cu ft 

p w mass density of water, 1.970 slugs/cu ft 

Subscripts: 
a air 

c capture station 

d diffuser exit 

e exit station 

f frontal 

i inlet 

m mixing chamber 

w water 

wa mixture of water and air 


L 

1 

2 

4 

9 


00 


free stream 



5 


I. EXPERIMENTAL INVESTIGATION OF HYDROKJCT 
By Elmo J. Mottard 

HYDRODUCT MODEL AND TEST APPARATUS 


The hydroduct model and towing apparatus are shown in figures 1 
and 2. The model was supplied by the manufacturer. The diffuser and 
nozzle sections were made of clear plastic to facilitate flow observa- 
tion. These two pieces were screwed onto a brass midsection which was 
welded to a strut. The strut had a duct which connected a compressed- 
air supply to the interior of the brass midsection. Metered air was 
injected into the stream of water through a cylindrical screen with 
openings of 0.001 inch in diameter. The space between the screen and 
the outer shell of the midsection constituted an annular settling 
chamber . 

Pressure-measurement orifices were provided at the maximum- diameter 
section of the diffuser and at three points around the inlet, located 
so that an axial asymmetry due to yaw, angle of attack, or distortion 
of the water surface could be detected. Provisions were made for keeping 
all lines to pressure taps completely filled with either air or water. 
Ityriamometer measurements were made of the axial force on the model. Nose 
and tail fairings shaped as shown in figure 1 were used for tare meas- 
urements. The towing tank in which the tests were made is described in 
reference 3* 


PROCEDURE 


The speed and airflow rate were held constant during the tests. 

In order to minimize strait tares, the tests were made at a depth of 
only 1 inch to the top of the hydroduct body , with the exception of a 
few tests at a 7-inch depth to investigate surface effects. Speeds 
of 20, 30 , 40, 50, and 60 feet per second and airflow rates of approxi- 
mately 0.1, 0.2, and 0.3 pound per second were used. A reference drag of 
the strut and hydroduct was obtained with the nose and tail fairings in 
place. 


RESULTS AND DISCUSSION 


The axial-force coefficient is plotted as a function of the ratio 
of weight airflow to free-stream velocity in figure 3. It is evident 
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that the hydroduct is capable of producing thrust since positive values 
of the axial-force coefficient were obtained. The increase in draft 
from 1 inch to 7 inches produced a constant decrement in positive axial- 
force coefficient of 0.07, which is attributed to the increased strut 
drag. It is believed, therefore, that no surface effect, such as a 
reduction in ram pressure or change in flow at the exit due to proximity 
of the free surface, occurred. 

The thrust coefficient, obtained by subtracting the drag coefficient 
of the body with the openings faired over from the axial-force coefficient, 
is presented in figure 4 . It is apparent from figures J and 4 that the 
force coefficients are functions primarily cf the ratio of weight air- 
flow to free-stream velocity. 

The variation of the pressure coefficient at the inlet with the 
ratio of weight airflow to free-stream velocity is presented in figure 5. 
The variation of pressure around the inlet was small, an indication 
that the flow was nearly symmetrical radially despite the proximity of 
the free surface. The pressure coefficient at the diffuser exit is pre- 
sented in figure 6. In a manner similar to the force coefficients, the 
pressure coefficients (figs. 5 and 6) are functions primarily of the ratio 
of weight airflow to free-stream velocity. The inlet velocity ratios 
were computed from the measured inlet pressures and are presented in 
figure 7 , That the maximum inlet velocity ratio without airflow (1.53) 
approached the ratio of exit area to inlet area (I.65) indicates that 
the system total-pressure losses were small for these conditions. For 
free-stream velocities greater than 40 feet per second, cavitation 
occurred at the inner lip sections of the inlet which precluded deter- 
mination of the correct inlet velocity ratio from inlet static-pressure 
measurements. The inlet velocity decreases with increasing airflow 
until air starts to "puff" out of the inlet; this occurs with a ratio 
of weight airflow to free-stream velocity greater than about 0.008. 

The inlet velocity can no longer be determined when air leaks out the 
inlet . 

In figure 8 is presented the variation of propulsive efficiency 
with thrust coefficient. The propulsive efficiency is defined by 


Efficiency - 


Forward speed x Eet jet thrust 
Power available from expanding air 


where the denominator is the power ideally available from an isothermal 
expansion of the existing airflow from the pressure measured in the 
diffuser to the free-stream pressure. An isothermal expansion was 
selected as a criterion approachable by the actual expansion because 


H CVI -=t CT\ 



J H OJ4- Q\ 
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of the great heat capacity of the water (relative to that of air) and 
the rapid heat-transfer rate attainable due to the mixing. Use of an 
adiabatic expansion as a criterion, for example, could result in effi- 
ciencies exceeding 100 percent. The maximum efficiency obtained in 
these tests was approximately 53 percent at a thrust coefficient of 0.10 
and speed of 60 feet per second. The range of test variables, however, 
was insufficient to determine the maximum attainable efficiency as a 
function of thrust coefficient or free-stream velocity. The data indi- 
cated that the maximum efficiency probably occurs at a free-stream veloc- 
ity greater than 60 feet per second (fig. 8). The hydroduct is compar- 
able to a supercavitating propeller for application to a high-speed 
hydrofoil configuration; however, the thrust coefficient of 0.10 was 
small compared with that of 1.07 for a supercavitating propeller (com- 
puted from ref. k ) . 

Consideration has been given to the possibility that exit choking 
has an effect on thrust production and efficiency. Examination of ref- 
erence 5 indicated that air-water mixtures of the range of the investi- 
gation could have sonic velocities as low as 70 feet per second. In 
the absence of exit static-pressure measurements, exit velocities were 
calculated by assuming that the exit pressure was equal to the free- 
stream static pressure; most of the resultant values were equal to or 
greater than 70 feet per second, indicating the possibility of exit 
choking. If the exit flow was indeed choked, the exit pressures might 
be greater than the ambient pressure so that some additional expansion 
of the flow through the use of a larger exit area may produce greater 
thrust . 


II. THEORETICAL INVESTIGATION OF AN 
AIR-POWERED UNDERWATER RAMJET 
By Charles J. Shoemaker 


ANALYSIS 


This theoretical development is for an underwater ramjet that uses 
air as a "fuel. 11 The ramjet cycle is as follows: the water captured by 

the inlet passes through a diffuser into a mixing chamber; here, com- 
pressed air is injected at the mixing-chamber pressure and water tem- 
perature; the resulting water-air mixture is discharged through an exit 
nozzle. Because the water-air mixture is less dense than pure water, 
the exit velocity exceeds the forward speed and a thrust results. 
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The procedure to be used in the theoretical development of an 
underwater ramjet cycle is as follows: A basic equation for thrust 

coefficient is presented and this basic equation is then evaluated by 
use of simplifying assumptions concerning the fluid properties. Two 
assumptions are made concerning the type of flow through the exit nozzle: 
(l) that the flow is incompressible and (2) that the flow is compressible . 


Basic Thrust-Coefficient Equation 

A simple underwater ramjet configuration is shown in sketch 1, as 
follows : 



Sketch 1 


The internal thrust of a ramjet cycle for steady flow is determined 
by the rate of change of momentum of the fluid passing through the ramjet. 
If the free-stream pressures of the air and water are assumed to be equal, 
then the internal thrust of the configuration shown in sketch 1 is 

T = “w( v e “ Voo] + ma(v e - vj) + (p e - p m ^A e 


where m^ and 1 % are, respectively, the mass flow per unit time of 

water and of air entering the ramjet. By rearranging terms and using 
the equality = p w VooA c , the internal thrust becomes 


T 




^ H OJ 4" 



^ H OJ -=t Q\ 
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The thrust coefficient based on capture area is defined as 


I P W V oo Ac 

Upon substituting the thrust relation into the preceding equation and 
canceling terms , the following basic equation for internal thrust coef- 
ficient is obtained: 


c T,c = 2 



( 1 ) 


Subsonic Discharge 

A basic assumption of this analysis is that the ramjet operates 
with a subsonic exit. For a subsonic exit, it can be demonstrated that 
the exit and free -stream pressures are equal - that is, p e = p^. Then 
examination of equation (l) shows that the term containing the pressure 
and areas becomes zero. 

To evaluate the thrust coefficient requires the determination of 
the velocity ratio V e /Vco. This ratio is most easily evaluated by 

studying, in succession, the processes of induction, mixing, and dis- 
charge. The assumptions for the induction and mixing processes remain 
the same regardless of whether the exit nozzle flow is treated as incom- 
pressible or compressible. 

Inlet and diffuser .- In the induction process it is assumed that 
the water acts as an incompressible inviscid fluid. Also, the diffuser 
is assumed to deliver the water to the mixing chamber at a negligible 
velocity. With the use of these assumptions and Bernoulli's equation, 
the total pressure in the mixing chamber can be expressed as 


P t,m = P - + | P v V °° 2 


( 2 ) 


Mixing of air with water .- The assumptions of the mixing process 
are: the air is injected at the total pressure and total temperature 

of the water, the air and water mix homogeneously without loss in total 
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pressure, the air does not dissolve in the waler nor does the water form 
a vapor, and the air and water components always have the same velocity. 

Incompressible flow in exit nozzle .- One method of analyzing the 
flow through the exit nozzle is to treat the mixture of air and water 
as an incompressible fluid. The nozzle total pressure is assumed to be 
a constant and to be equal to the mixing- chamber pressure. 

Velocity in relation to density: Bernoulli ' s equation for incom- 

pressible flow evaluated at the nozzle exit is 


1 .. 2 

Pt,m - Pe + 2 p wa,e e 

Substituting the mixing-chamber total pressure; (eq. (2)) into this equa- 
tion and using the relation = p e results in 


Poo + 2 Pw^ 00 Poo + 2 p wa,e^e 
Solving now for the velocity ratio V e /Voo gives 



(3) 


Density ratio: The following development gives a general expres- 

sion for the density ratio. The volume flow of the water-air mixture 
at a given station is 

Volume flow = ^ = K + m a) = + 'ZsS 

Pwa Pwa \Pw p a J 


Solving for l/p wa gives 


^ H CM G\ 



H CVJ -=f ON 
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tot + M jJ_l + Ea 

1 _ \Py £aj _. ^yyPw ^w p a/ 

p wa (% + %) _1 A + »a\ 

m w\ %/ 


If both sides of the equation are multiplied by p w and then simplified, 
the following general expression is obtained for the density ratio: 




+ ^a M 

rcy Pa/ 



(M 


Note that the air and the mixture densities are associated values and 
must be evaluated at a common point. 

Velocity ratio: The complete expression for the velocity ratio in 

an incompressible mixture of air and water is obtained from equations (3) 
and ( 4 ) and is 



Conqpressible flow in exit nozzle .- Another method of analyzing the 
flow through the exit nozzle is to treat the water-air mixture as a com- 
pressible fluid. This method is preferred because the actual mixture 
is compressible. 

For this development, the nozzle total pressure is again assumed 
to be a constant and equal to the mixing-chamber pressure and the mixture 
of air and water is assumed to act as a compressible gas. By this 
approach, the velocity ratio is obtained as a function of an appropriate 
adiabatic exponent 7. The estimation of 7 is discussed in detail in 
a subsequent section. 

Velocity in relation to density: In order to obtain the velocity 

ratio V e /Voo, the following compressible-flow relation is used: 
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^5= (l^Me 2 ) 

Subtracting 1 from each side of this equation and inserting 

M e 2 = V e 2 /a e 2 = V e 2 /( 7Pe/pva,e) S ives 



with the assumption that Pt m = Pt,wa> ^' Qe mixing- chamber total pres- 
sure obtained from equation (2) can be rearranged to give 

%o _ P t ,wa _ Using this relation in the preceding equation and 
Poo Poo 

solving for V e produces 



Dividing both sides of this equation by V M and inserting p w /p w under 
the radical gives 



1 2 

When p e = p m is substituted and the identity q^ = ^ P w v °° is used 

in the preceding equation, rearranging the resulting equation produces 
the following expression for the velocity ratio: 


^ H C\l -3* ON 



J rH OJ4- On 
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This velocity ratio (eq. (6)) is the incompressible-! low relation 
(eq. (3)) multiplied by a function to compensate for compressibility. 

The term under the radical which is enclosed by the braces is referred 
to herein as a compressibility function and is designated f(7) because 
it contains an unknown 7. 


Velocity ratio: The substitution of equation ( 4 ) into equation (6) 

gives the complete velocity ratio for compressible flow, which is 



( 7 ) 


Estimation of adiabatic exponent: To obtain the adiabatic 

exponent 7 needed in equation (7)j the following relation is used: 


-( iQ g e p) 
■ (l°g e p) 


= 7 


( 8 ) 


Equation (8) can be obtained from the speed-of- sound equation or from 
the relation p = Constant x p 7 . The latter relation is used in this 

development. Differentiation of this relation yields 

/p7 \ 

dp = Constant X 7I— jdp 


Substituting p for the quantity Constant x p 7 and rearranging gives 



dp 

P 



Ik 


Using the relations d(log e p) 


dp 

P 


and 


d(lcg e p) = y- 


produces 


d(log e p) 

— = y 

d(log e p) 


The evaluation of 7 is begun by differentiating equation (4) 
with p wa and p a as the variables. When equation (4) is inverted 

and differentiated, the result is 


dp. 


wa 


w 


f %\f%Pw^Pa] 
y + m-^a 2 )_ 

(1 + Hi 

\ %P a 


Multiply both sides by P w / P wa to obtain 


dp. 


wa 


w L + ^VMfPwVfPa 

. m w/\ \ p a'\ p a 


Pwa 


wa 


l HV P a j 


Replacing Py-^Pwa "by equation (4) and canceling like terms gives 


^Pwa 

Pwa 


m a Pw 
% Pa 
+ Pw 

% P^ 



However, 


^ = d( loge p) 


and therefore 


d(l°g e Pwa) = d ( lo 6e P a ) 


/ 

^w Pa 




m a Pw 
1 + - 

m w Pa/ 


vo -r ro m r* 



^ H OJ ON 
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Invert this relation and multiply by d(log e p) to obtain 


d(log e p) = d( log e p) 
d(l°g e P wa ) d ( lQ ge Pa) 


By use of equation (8) and rearranging terms, the following general 
relation for 7 is found: 


,1 + £si 

m w Pa 


m a Pw 
Pa 


7 


^w Pa 
m a Pw 


(9) 


As the mixture travels through the nozzle , 7 is continually 

changing. Therefore, the 7 used to evaluate equation (7) was an aver- 
age of the mixing-chamber and exit values. The average value of 7 is 


7 


7 wa,e + 7 m 
2 


Substituting for 7 wa ^ e and ?m using equation (9) and rearranging 
terms produces 


7 


r“ 


^ m 3. Pw y 



If it is assumed that the air expands isothermally , the relation 
p = p^RT shows that p is proportional to p. Then, noting that 

Pm = Pt = Poo + q<» y ields 


7 = 7 P 


1 + 1 ^w Pa^e L + + 

2 ma p w \ p e / 
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For the assumption that the nozzle is not choked, the free-stream pres- 
sure equals the exit pressure (p^ = p e ). Making this substitution in 
the preceding equation and combining terms gi’ r es the expression for 
the average adiabatic exponent: 


7 


+ ^ Pa > e (i + i ^°\ 
% P w ^ 2 P oo) 


( 10 ) 


In the differentiation leading to equations (?) and (10), p w was 
assumed to be constant. This is a reasonable assumption for the thrust 
calculations. However, if it is desired to have an expression for the 
speed of sound which covers the range from air only to water only, then 
the compressibility of water must be included. 

Component efficiencies .- Throughout this theoretical development, 
the maximum thrust coefficient was obtained by assuming that the component 
losses were negligible. In this section, a procedure for including 
the effect of the various component losses is given. When the analytical 
and measured thrust coefficients are compared, the losses can be 
estimated. 

The effect of the component losses can be expressed as efficiencies 
in the following manner. Let the inlet efficiency be defined so 

that the mixing- chamber total pressure is 


Pt,m = Poo + p w v« 2 

Using this inlet total pressure and solving for the velocity ratio, by 
the method used in the section entitled "Incompressible flow in exit 
nozzle," gives a relation similar to equation (3)- The velocity ratio 
is now 



The combined efficiencies of the inlet, mixing process, and exit nozzle 
are defined as an overall efficiency g(r|). With g(rj), the incom- 
pressible velocity ratio is then 


vo ro h tr* 



VO 4 =- ro H f 
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By analogy, equation (6) which is for compressible nozzle flow becomes 



( 11 ) 


Working Form of Thrust-Coefficient Equation 

The assumptions made to obtain the working form of the thrust- 
coefficient equation which was used to study the experimental data are 
again stated. 

Throughout the ramjet cycle, the component losses were neglected. 

The water was considered to be an incompressible fluid. In the mixing 
chamber the water velocity was negligible. The injection of air 
resulted in a uniform mixture that did not separate. This mixture was 
discharged through a nozzle with a subsonic velocity. The nozzle total 
pressure was constant and equal to the mixing-chamber pressure. Static 
pressure at the nozzle exit was equal to the free-stream static pressure. 
The mi xture behaved as a compressible gas. An evaluation of the average 
synthetic 7, based on an isothermal expansion, was made for this gas. 


flow 


C T,c 


The internal thrust coefficient for subsonic compressible nozzle 
obtained from equations (l) and (7) is 



An average value of 


7 } given by equation (10), was used. 
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Sonic Discharge Considerations 

As mentioned previously for subsonic discharge, the nozzle exit 
pressure and free-stream pressure are equal. If the exit velocity is 
equal to sonic velocity, the exit and free-stream pressures are not 
necessarily equal and a different procedure is necessary to evaluate 
the thrust coefficient. Herein are given two ways of estimating the 
sonic velocity of a water-air mixture. Then, a suggested procedure to 
determine the thrust for a sonic discharge is presented. 

Estimation of sonic velocity .- The speed of sound of the mixture L 

is a function of the component densities which are determined by their 1 

respective temperatures. There are two expansion processes which 
bracket the temperature range. One is an isentropic expansion and the 
second is an isothermal expansion. The speed of sound obtained from 
each of these expansions is presented and compared. 

Isentropic expansion: If the water-air mixture is assumed to pass 

through the nozzle so quickly that there is kittle time for heat trans- 
fer, the expansion process is essentially adiabatic. Also, the losses 
in a nozzle are assumed to be negligible since they are usually small. 

With these assumptions, the expansion process is isentropic. 

A relation for the speed of sound is obtained by use of the fol- 
lowing equation for the velocity of sound: 


1 _ dp 

a 2 dp 


( 13 ) 


The right-hand side of equation (lj) is evaluated by the use of the 
mixture density obtained from equation ( b ) which is 


p wa 





% Pa J 


Here the water is to be treated as a compressible fluid; then p w , P wa ; 
and p a are associated variables and must be evaluated at a common 

point. When both sides of the equation are <Livided by a common pressure 
and differentiated with respect to this pressure, the following relation 
is obtained by rearranging terms: 


OJ U 03 
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L 

1 

2 

4 

9 


^wa 

dp 


+ 



+ ®a 

m w c>p 



Using equation (13) in the preceding equation yields 


^a 


- i + 
l “vy 


p w) a 2 
' a w 


1 + m a 1 


m. 


V Uq 


w 


^a\ 

m T . r 


Solving for a wa produces the following general relation for the speed 
of sound in the mixture: 



For mass-flow ratios of the data which were sufficient to produce 
thrust, the denominator of equation (14) is essentially unity. 

If the water had been assumed to be incompressible, the speed-of- 
sound equation would be 


a wa 


a a 


m,. 


+ nig. 


1 + 


v“aAPW 


(15) 


Equation ( 15 ) can be reduced to the numerator of equation (14). 

Isothermal expansion: Another method to obtain the speed-of-sound 

equation is given by Prandtl in reference 5- Prandtl assumed the air 
temperature to be equal to the practically constant water temperature. 
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The water density was assumed to be constant. Also, the total mass flow 
was chosen to be unity. The speed-of -sound relation given by Prandtl 
is equivalent to 



(16) 


Comparison of speed-of-sound relations: In order to compare the 

two speed-of-sound relations , equation ( 15 ) was modified. Substituting 

a a = \Jva(p/pa) anh % + = 1 into equation (15) and rearranging 

gives 


a wa 



(17) 


When water is assumed to be incompressible, the sonic velocity of 
a water-air mixture for an isothermal expansion (eq. (16)) is less 
than the sonic velocity for an isentropic expansion (eq. (17)) by the 

factor \Jy 

Choked flow .- When sonic velocity is obtained at the exit of a 
converging nozzle, there are two distinguishing features about the flow: 
First, the pressure at the exit is generally unequal to the free-stream 
pressure and, second, the nozzle is discharging its maximum mass flow 
(choked flow). This choking of the exit restricts the total mass passing 
through the ramjet. 

In the analysis of a ramjet cycle with choked exit flow, the assump- 
tions for the inlet, diffusion, and mixing processes are the same as 
those for nonchoked exit flow. The mixing-chamber velocity is again 
assumed to be negligible. If the mixture of water and air is assumed 
to behave as a compressible gas, the mixture properties can be evaluated 
as a function of an adiabatic exponent 7. By assuming an exit pressure 
p e that is equal to or greater than p^, 7 can be obtained from equa- 

tion ( 9 )* For an exit Mach number of unity and with 7, the other var- 
iables in equation (l) can be calculated by using the compressible gas 
relation for P+/p, Bernoulli* s equation, and the continuity equation. 


vo -p- ro t— ' t"* 



J H OJ ON 
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In a comparison of the nozzle exit velocity (eq. (7)) and th e 
sonic velocity of the water-air mixture (eq. (l4)) , it was found that 
all of the forward-velocity data at 50 and 60 feet per second had 
reached or slightly exceeded sonic velocity at the nozzle exit. Inas- 
much as the theoretical thrust coefficients computed by equation ( 12 ) 
are not for choked flow, a separate analysis is required for these data. 
However, an exploratory computation for a forward velocity of 60 feet 
per second and a mass-flow ratio of 0.02 indicates that the thrust coef- 
ficient for choked flow is only about 3 percent more than that for the 
subsonic flow. 


Determination of Thrust Coefficient 

The computations for the theoretical thrust coefficient are based 
on compressible nozzle flow with a subsonic exit velocity. 

Compressibility function .- In order to evaluate the compressibility 
function f ( 7 ) used in equation ( 6 ), the density of the air at the exit 
is needed. This density was obtained by assuming that the air temperature 
was equal to the free-stream water temperature of 505 ° R. 

Figure 9 is a plot of f(z) as a function of q^/Poo for several 
values of 1 %^%. The compressibility functions calculated for the 
particular values of 1 %/my obtained from the data of part I are 
included in the figure. 

Theoretical internal thrust coefficient .- Because the experimental 
thrust coefficients are based on frontal area, the theoretical thrust 
coefficients have been based on frontal area by multiplying equation ( 12 ) 
by the area ratios ( A c /A e ) and (A e /Af). The first area ratio is obtained 
from the continuity equation and is 


Ac 

A-e 


The second area ratio Ae/Af is obtained from the model dimensions. 
Figure 10 shows the relation between the internal thrust coefficient 
based on frontal area and the mass -flow ratio 1 %/my- for several values 
of f( 7 ). The thrust coefficients evaluated for the particular mass- 
flow ratios obtained from the data of part I are plotted in this fig- 
ure, with the values of f( 7 ) indicated beside them. 


Ve/Vo 
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Figure 9 can be used in conjunction with figure 10 to provide an 
estimation of the theoretical thrust coefficient. By choosing a for- 
ward velocity of the ramjet , the ratio can be calculated. Then, 

for an assumed mass-flow ratio figure 9 gives the theoretical 

compressibility function f(7). With f(7) and the assumed m^jm^y 

from figure 10 the corresponding theoretical thrust coefficient can be 
estimated. 


COMPARISON OF THEORY WITH EXPERIMENT 


Thrust Coefficient 

The theoretical and measured thrust coefficients are presented in 
figure 11. The theoretical values have been computed by using the mass 
flow obtained from the data of part I and have been faired according to 
values of constant forward velocity. The differences in the measured 
and theoretical thrust coefficients are due to three factors: (l) The 

experimental method of obtaining the measured thrust coefficient does 
not give the exact internal thrust coefficient. (2) There are internal 
losses associated with the inlet, the mixing process, and the exhaust 
nozzle. (3) The theory depends on the validity of the assumptions used 
to describe the mixture properties. With the viewpoint that the effect 
of factor (l) is negligible and the assumptions of factor (3) are valid, 
the internal losses become responsible for the differences in the theo- 
retical and measured thrust coefficients. 


Efficiency 


When the efficiencies of the ramjet components are taken into _ 
consideration , the curves of figure 10 represent the product j~f( 7 )J j^g(r])J 


as indicated by equation (ll). The value of this product is obtained 
by transferring the data points of figure 11 to figure 10. It is seen 

that most of the data lie close to a value of [f(7)][g(l)] = 0.6. An 

estimation of the overall efficiency is obtained by dividing the pre- 
vious relation by the jf( 7 )j of the data whi:h is given in figure 10. 
Then, by using an attainable inlet efficiency, the mixing and exit 
nozzle efficiencies can be estimated. 


rH c\j -=t o\ 
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In the estimation of the inlet efficiency; it is assumed there is 
isentropic compression between the capture and inlet stations. For the 
remaining compression, regular subsonic diffuser losses are used. From 
subsonic diffuser data, these regular subsonic diffuser losses are 
known to depend on the diffuser wall angle and to be a certain percent 
of the losses for a sudden expansion to the same area ratio. For the 
inlet configuration used in the experiment, the efficiency is estimated 
to be 





The solid curves of figure 12 give the overall efficiency as a 
function of The dashed curves in figure 12 are the results of 

dividing the overall efficiency by the estimated inlet efficiency; this 
gives an estimation of the mixing and nozzle efficiency. Inasmuch as 
the nozzle efficiency is usually very close to unity, these dashed 
curves could be assumed to be the mixing efficiency. 


CONCLUDING REMARKS 


It has been found from a preliminary experimental investigation of 
an underwater ramjet that a positive thrust-minus-drag force was pro- 
duced at every test speed. The force and pressure coefficients were 
functions primarily of the ratio of weight airflow to free-stream veloc 
ity for the conditions of the investigation. The maximum propulsive 
efficiency based on the net internal thrust and an isothermal expansion 
of the air was approximately 53 percent at a thrust coefficient of 0.10 
The performance of the test model may have been influenced by choking 
of the exit flow. A more complete study is required to determine the 
applicability of the hydroduct. 

The basic assumption of the theoretical analysis is that a mixture 
of air and water can be treated as a compressible gas. However, more 
information on the properties of air-water mixtures is required to 
confirm this assumption or to suggest another approach. A method is 
suggested from which a more complete theoretical development, with the 
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effects of choking included, may he obtained. An exploratory computa- 
tion, in which this suggested method was used, indicated that the 
effect of choked flow on the thrust coefficient was minor. 


Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Air Force Base, Va., September 28, 1961. 
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Towing carriage 



Figure 1. - Schematic sketch of hydroduct and test setup. Linear dimensions are in inches 
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Figure 2.- Photograph of hydroduct model. L-60-1921.1 



Air blowing out of inlet 



Figure 3*- Variation of the axial- force coefficient with the ratio of 
the weight airflow to the free- stream velocity. 
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Figure 4.- Variation of the thrust coefficient based on frontal area 
with the ratio of weight airflow to f ree-stream velocity. 





Pressure coefficient at the inlet, 



Figure 5*- Variation of the pressure coefficient at the inlet with the 
ratio of weight airflow to free-stream velocity. 


Figure 7 .- Variation of the inlet velocity ratio with the ratio of 
weight airflow to free-stream velocity. 



Efficiency (based on an ideal isothermal expansion) 
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Figure 9.- Compressibility function. m a based on mixing- chamber pressure and 505' 
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Figure 10.- Theoretical thrust coefficient based on frontal area, based on mixing- chamber 

static pressure and 505° R* 
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Figure 11.- Comparison of theoretical and measured thrust coefficients. 
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Figure 12, - Overall efficiency and mixing efficiency. 



